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SI Experimental methods 

This section gives additional details about our strains and experimental methods. 
Sl.l Strains 

The strains used in this work are listed in Table SI. The mutualistic strains Leu FBR Trp~ and Leu - Trp FBR 
are W303 strains that are prototrophic for all amino acids except for either leucine or tryptophan. They all 
have the genetic background W303 MAT a canl-100 hmlaA::BLE leu9A::KANMX6 prACTl-yCerulean- 
tADHl@URA3, and differ only at the LEU4 and TRP2 loci, as well as in their fluorescent marker at the 
HISS locus (see below). The strains contain the S288C version of BUD4 (standard W303 strains have a 
mutated BUD4 [1], which causes alteration in the budding pattern and occasional failures in cell separation 
after cytokinesis). The deletion of the silent mating type locus HMLa prevents mating-type switching, and 
therefore mating, since all our strains are MATsl at the mating type locus. The two strains thus behave 
essentially as two different species. Leu9 and Leu4 are isozymes for the first step in leucine biosynthesis [2], 
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Strain name 


Genotype at LEU 4 and TRP 2 loci 


fluorescent proteins 


Leu FBR Trp - yMM60 
Leu" Trp FBR yMM65 


LEU4 FBR trp2A::NATMX4 
leu4A::HPHMX4 TRP2 FBR 


ymCitrine, yCerulean 
ymCherry, yCerulean 


yMM61 
yMMo4 


LEU4 Fati trp2A::NATMX4 
leu4i-^.::nPHMA4 1 RP2 


ymCherry, yCerulean 
ymCitrine, yCerulean 


yMM26 
yiVllvlal 
yMM29 
yMM32 


LEU4 hali tr P 2A::NATMX4 

LEU 4 trp2A::NATMX4 
leu4A::HPHMX4 TRP2 


ymCitrine 
ymCherry 
ymCherry 
ymCitrine 


non-interacting yJHKlll 
strains yJHK112 


LEU 4 TRP 2 
LEU 4 TRP 2 


ymCitrine 
ymCherry 



Table SI: Strains used in this work. All strains have the genetic background W303 MATsl canl-100. All 

yMM strains have in addition hmlaAr.BLE leu9A::KANMX6. 

The fluorescent markers are incorporated as follows: 

ymCitrine: hisSA ::prA CTl-ymCitrine-tADHl :HIS3MX6, 

ymCherry: hisSA ::prA CTl-ymCherry-tADHl :HIS3MX6, 

yCerulean: prACTl-yCerulean-tADHl@URA3. 

which makes the leu9 deletion necessary in order to make a leu4 -deletion auxotrophic for leucine. With 
the mCerulean gene under the actin promoter prACTl, all strains constitutively express a cyan fluorescent 
protein. We checked by flow cytometry that the cyan fluorescence intensity profiles were indistinguishable 
for all strains. 

The leucine overproducing strain Lcu FBR Trp~ is the strain yMM60, which has the (additional) genetic 
modifications his3A::prACTl-ymCitrme-tADHl:HIS3MX6 LEU4 FBR trp 2 A : :NA TMX4 . It constitutively 
expresses the yellow fluorescent protein mCitrine and is a tryptophan-auxotroph. The LEUA FBR allele 
differs from the wild-type LEU4 by the deletion of the codon 548; it is functional in leucine biosynthe- 
sis, but feedback-resistant (FBR) to inhibition by the end-product leucine. It was identified in Ref. [3] 
(named there LEU4-1)- This strain is depicted as yellow (in accordance with its fluorescent protein) in the 
main text. The tryptophan overproducing strain Leu - Trp FBR is the strain yMM65, with his3A::prACTl- 
ymCherry-tADHl:HIS3MX6 leu4A::HPHMX4 TRP2 FBR . It constitutively expresses the red fluorescent 
protein mCherry and is a leucine-auxotroph. The TRP2 FBR gene is the allele TRP2 — S76L, which was 
identified in Ref. [4] (designated L^q in this reference) as a feedback- resistant, functional version of TRP2. 
This strain is depicted as blue in the main text (to enhance its contrast with the yellow color of the partner 
strain) . 

To rule out effects of the fluorescent proteins on the fitness or other phenotypes of the 
strains, we also used strains with 'swapped' fluorescent proteins: yMM61, with his3A::prACTl- 
ymCherry-tADHl:HIS3MX6 LEU4 FBR trp2A : :NA TMX4 , and yMM64, with his3A::prACTl-ymCitrine- 
tADHl:HIS3MX6 leu4A::HPHMX4 TRP2 FBR . In order to compare feedback-resistant and wild-type amino 
acid production we used the following strains: the FBR strains yMM26 and yMM31, which have the same 
genetic background as yMM60 and yMM65, respectively, except for being URA3 wild-type at the Ura3 locus; 
and the WT strains yMM29 and yMM32, which are LEU 4 his3A::prACTl-ymCherry-tADHl:HIS3MX6 and 
TRP2 his3A::prACTl-ymCitmne-tADHl:HIS3MX6, respectively. 

SI. 2 Fitness costs of amino acid overproduction due to feedback resistance 

The relative fitness in direct competition in liquid culture was measured by a flow-cytometer-based compe- 
tition assay as described in Ref. [5] . The relative fitness on plates was determined by analyzing the shape of 
sectors of the two competing strains in a colony as described in Ref. [6] . 

We first checked that the fluorescent proteins mCitrine and mCherry were neutral with respect to each 
other both in the liquid and the plate competition assays by verifying that the strains yMM60 and yMM61 
as well as yMM63 and yMM64, which differ only in their fluorescent protein, have the same fitness in CSM 
(for medium definitions, see the Materials and Methods section of the main text). 

We competed the mutualistic strains Leu FBR Trp~ (yMM60) and Lcu~ Trp FBR (yMM65) in liquid CSM, 
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Figure SI: Colony expansion velocities. A) Radius increase over time for individual colonics of the mutu- 
alistic strains Leu - Trp FBR and Leu FBR Trp~, inoculated at different start fractions /q of blue cells, under 
conditions of no (red dots) and obligate mutualism (green dots). In the no mutualism case, all colonies 
exhibit similar radius increases and reach a velocity of 20.7±0.4^m/h for times larger than 90h (upper red 
line). In contrast, obligately mutualistic colonies start expanding at different points in time: Colonies with 
start fraction /o = 0.15 start expanding at later times (solid purple line through green dots), and colonies 
with /o = 0.01 do not expand at all (solid red line through green dots). Expanding colonies reach a velocity 
of 10.7±0.6/im/h for times larger than 90h (upper green line), indicating that the mutualistic strains have 
reached steady-state expansion. B) Growth velocities for times larger than 90h are independent of the start 
fraction / 0 both for non-mutualistic (red circles) and for obligately mutualistic (green circles) colonies. Solid 
lines indicate the average velocities, which are the slopes of the corresponding straight upper red and green 
lines in (A). 



and found that Leu~ Trp FBR had a 2% fitness advantage compared to Leu FBR Trp~. The same was true in 
a competition of the strains yMM61 and yMM64 with 'swapped' fluorescent proteins. This fitness difference 
is entirely due to leucine feedback-resistance, as the Trp FBR strain yMM31 and the Trp WT strain yMM32 
have the same fitness, while the Leu FBR strain yMM26 and the Lcu WT strain yMM29 also have a 2% fitness 
difference. Presumably, leucine overproduction due to feedback-resistance is more costly than tryptophan 
overproduction because leucine is more abundant than tryptophan, with ~10% of yeast amino acids being 
leucine and only ~1% being tryptophan [7-9]. 

SI. 3 Radial expansion velocities 

Radial expansion velocities. When a yeast colony expands on an agar surface, its radius increases linearly 
with time after an initial transient of about 4 days, see Fig. S1A. In the linear regime, we characterize colony 
growth by an expansion velocity, defined as the slope of the radius versus time. For both mutualistic and 
non-mutualistic colonies, the steady-state expansion velocity is independent of its initial conditions, such as 
the inoculation volume and density or the initial mixing ratio of the blue and yellow strains. 

Initial transient in radius increase. In contrast, the transient towards the linear regime can depend 
on the initial condition, since it takes longer to reach steady-state growth mode when the initial conditions 
are further away from it. As shown in Figure 3 in the main text, obligately mutualistic colonies are driven 
towards an optimal fraction /* = 0.5 of blue cells. Thus, when colonies are inoculated with a start fraction 
fo close to /*, they approach the steady-state growth velocity faster than when inoculated with a start 
fraction further from /*. This effect is particularly noticeable for small fo, see the orange and red data in 
Fig. SlA, i.e. when there is a very small fraction of leucine requiring cells. In the case of fo = 0.01, the 
colony never grew on experimental timcscalcs. 

Growth as function of the amino acid concentrations. Figure 2A in the main text shows that, as 
a function of a required amino acid, the growth velocity increases linearly before it saturates at a plateau. 
This data is reproduced in Fig. S2 (blue diamonds in (A), yellow squares in (B)). In addition, this figure 
shows that the growth velocities do not depend on the concentration of the non-required amino acid (yellow 
squares in (A), blue diamonds in (B)). The maximal growth velocities (plateau values) for Leu - Trp FBR and 
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Figure S2: Comparison of growth velocities on plates and growth rates in liquid of Leu FBR Trp - (yellow 
squares for velocities and green circles for rates) and Leu - Trp FBR (blue diamonds for plates and red stars 
for liquid) for (A) varying leucine and (B) varying tryptophan concentrations. Growth velocities and rates 
are independent of non-required amino acids, but increase linearly with a required amino acid concentration 
before they saturate at a plateau. When the liquid growth rates (red stars for Leu - Trp FBR and green circles 
for Lcu FBR Trp - ) are scaled by a factor of b vg = 40 /im, they are very similar to the growth velocities on 
plates (blue diamonds for Leu - Trp FBR and yellow squares for Leu FBR Trp - ). However, crossovers from 
linear to plateau regimes for liquid growth occur at concentrations [leu] c = 762 fiM. (red vertical line) and 
[trp] c = 97.6 fjM (green vertical line) that are lower by a factor of 1 + K = 1.18 than the corresponding 
crossover concentrations (blue and yellow vertical lines) for the plate growth velocities, as explained in 
Section S2.2. 



Leu FBR Trp - are «l = 20.3 /im/h and ut = 19.8 /xm/h, respectively. The 2.5% difference in the plateau 
growth rates agrees, within experimental error, with the difference in fitness between the two strains as 
measured by the competitive fitness assay described in Section SI. 2. 

We also measured the growth rate in well-mixed culture as function of the external amino acid concen- 
tration. As shown in Fig. S2, the growth rates behave very similarly to the growth velocities. In fact, when 
scaled by the same factor of b vg = 40 /im, the growth rate plateaus for both Leu - Trp FBR and Leu FBR Trp - 
overlap with the growth velocity plateaus. The growth rates in the linear regimes are similar but not iden- 
tical, since the growth rates in liquid increase slightly faster than the growth velocities on solid medium. 
This leads to higher crossover concentrations from the linear to the plateau regime for plate growth (blue 
and yellow vertical lines) than for liquid growth (red and green vertical lines). We normalize the leucine 
and tryptophan concentrations such that the crossover concentrations arc at «el = [leu]/ [leu] c = 1 and 
n ET = [trp]/[trp] c = 1 for liquid growth, with [leu] c = 762 fjM and [trp] c = 97.6 /iM, since the liquid con- 
centrations reflect the concentrations felt by the cells. The slower growth on plates is due to an amino acid 
depletion layer around the colony and leads to crossover concentrations which are larger by a factor of 1.18, 
see also Section S2.2. 

SI. 4 Determination of yeast colony growth parameters 

This subsection summarizes how we determined the model parameters listed in Table S2. 

Cellular diffusion constant D s and genetic diffusion constant D g . We determined the spatial 
diffusion constant D s (for cellular diffusion due to the jostling of cell division) and the genetic diffusion 
constant D g (which characterizes the strength of genetic drift) with the sectoring assay described in Ref. [10]. 
When performing this assay with the non-interacting equally fit strains yJHKlll and yJHK112, we obtained 
Ds/tq = 15/im 2 /h and D g jr g = 1.3/iin/h, where r g = 1.5 h is the yeast generation time. The value 
of D s for cellular diffusion lies in the expected range: In each generation, cells push by a diameter of 
a = 5 /im, so we expect the diffusion constant to be of the order of D s = a 2 = 25 /im 2 per generation, 
or D s /r g = 17/im 2 /h. The value of the genetic diffusion constant D g allows us to estimate the effective 
population density p c w 1/^g [10]. The resulting effective density p c ss 0.5//im of about 2-3 cells in a 
region of the size of a cell is rather low, especially when considering that yeast colonies pile up to a height 
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of ~ 100 cells. This result is similar to even lower effective population densities for the bacteria E. coli 
and P. aeruginosa observed in Ref. [10], and implies that the dynamics of the boundaries between different 
genotypes are dominated by the behavior of the cells at the extreme edge of the growing colonies. 

Size b of the active layer. To estimate the size b of the active layer, within which cell proliferation 
drives radial growth, we use the observation that the maximal colony growth velocity v is b vg = 40 pm times 
the maximal growth rate g in liquid culture (Section SI. 3). Within a generation time r g ~ l/g, the number 
of cells increases by AN = Ubc, where U is the colony circumference and c is the cellular surface carrying 
capacity. In the same time r ff , the colony increases in area by AA = Uvr g = Uv/g, which corresponds to 
AiV = cAA new cells. Thus, v/g = b. With v/g = b vg = 40 /mi, we obtain the reasonable estimate of 
b = 40 fim for the size of the active layer [11, 12]. 

Surface carrying capacity c. Colonies grown between 4 and 7 days (to obtain a variety of sizes, and to 
measure in the same time window used for all other experiments) were imaged with a Zeiss Lumar Stereoscope 
to determine the colony area by image analysis (circle fit with MATLAB). Directly after imaging, a small 
agar pad with the colony on it was cut out of the plate and suspended in a known amount of PBS. All cells 
were washed off the agar by strong vortexing. The cell concentration in the PBS was then determined by 
using a Beckman Coulter counter, from which the total number of cells in the colony could be calculated. 
The surface carrying capacity was determined as the slope of the cell number versus colony area, with the 
result of 10 cclls//xm 2 . Taking into account that yeast colonies have a height of 0.5 — 1mm, this surface 
density corresponds to a volume density of about 0.01 cells//xm 3 , which is similar to dense spherical packing 
of yeast cells of diameter 5 /im. Since we are interested in the cell density at the colony boundary, and since 
yeast colonies are higher in the middle than at the boundary, we take c = 5cells//xm 2 [11]. 

Amino acid diffusion constant D a . Diffusion through the agar medium is similar to diffusion in water, 
since the agar pore size of about 500A [13] is much larger than the amino acid hydrodynamic radius of about 
3A [14]. The diffusion constant for the similar-sized leucine and tryptophan molecules is D & = 3mm 2 /h [15]. 

Diffusive flux rate d. The derivation of d = 5/h is described in Section S5. 

Amino acid crossover concentrations [leu] c and [trp] c . The leucine and tryptophan crossover 
concentrations between the regime where growth depends linearly on the amino acid concentrations and the 
regime where growth is saturated are determined in Fig. S2 as [leu] c = 762 pM = 2.3 ■ 10 9 molccules//xm 2 
and [trp] c = 97.6 pM = 2.9 • 10 8 molcculcs/^tm 2 . Here, wc have used the agar height of 5 mm to transform 
the three-dimensional volume concentrations to two-dimensional surface concentrations, i.e. to the number 
of molecules available to cells from the agar directly beneath. 

Amino acid uptake rates. From the literature, we estimate the amino acid uptake rates as k^ = 
10 5 molecules/s/cell for leucine [8,16,17] and k^ = 1.5 • 10 4 molecules/s/cell for tryptophan [8,18]. These 
values are in agreement with the value of the parameter k = k^ c/[leu] c /d = fcx c/[trp] c /rf = 0.18 (determined 
independently from our growth rate and velocity measurements in Fig. S2 in Section S2.2) and our above 
estimates of [leu] c , [trp] c , c and d. 

Amino acid secretion rates. We expect the amino acid leakage rates to be proportional to the 
internal cellular amino acid concentrations. Since yeast cells contain an order of magnitude more leucine 
than tryptophan [7-9], we expect the leucine secretion rate to be an order of magnitude larger that the 
tryptophan secretion rate. From our experimental result that our yeast strains grow twice as fast under non- 
mutualistic than under obligate mutualistic conditions (Section S3.1), wc determined the parameter p = 1.09 
from the equation p = 1 + k/2. With p — tl c/[lcu] c /d = ?'t c/[trp] c /d, we obtain the leucine and tryptophan 
secretion rates = 7 • 10 5 molecules/s/cell and rj- = 9 • 10 4 molecules/s/cell. These values are of the order 
of magnitude of amino acid secretion rates measured in Ref. [19], but are on the larger side, presumably 
because our strains are overproducing leucine and tryptophan due to the engineered feedback-resistance and 
thus leak large amounts of excess amino acids. 

S2 Modeling growth of the mutualistic yeast strains in the pres- 
ence of amino acids 

In this section we describe a model for the growth of our mutualistic yeast strains in the presence of amino 
acids in more detail than in the main text. The model parameters arc listed in Table S2. 
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Symbol 


Definition 


Value 


Description 


K 

P 


c A:t c 

[leu] c d [trp] c d 

r"L c tt c 

[lcu] c d [trp] c d 


0.18 
1.09 


reduced amino acid uptake rate 
reduced amino acid secretion rate 


k 
lb 


1 D a /d 
V l+« 

2D a 
b 


710 /Ltm 
0.83 /j,m 


length scale of amino acid depletion 
length scale of patch boundary diffusion 


s c 


Dl/D s 


0.18 


critical selection strength to overcome genetic drift 



Table S2: Parameters used in the theoretical calculations of Section S3 and S4. As given by the definitions in 
the table, these 'reduced' parameters are compounds of directly biologically relevant parameters: the leucine 
and tryptophan uptake rates /cl and fcx, the leucine secretion rates tl and fx, the leucine and tryptophan 
crossover concentrations [lcu] c and [trp] c , the cellular surface carrying capacity c, the amino acid diffusive 
flux rate d, the amino acid diffusion constant D a , the size b of the actively growing layer, the cellular diffusion 
constant D s and the genetic diffusion constant D g . The determination of the parameter values is described 
in Section SI. 4. 

S2.1 Model for amino acid dynamics 

Growth velocities. As shown in Section SI. 3, the growth rates and velocities increase roughly linearly 
with the concentration of the required amino acid until they saturate at a plateau. This motivates us 
to write the growth velocities as function of the rescaled amino acid concentrations til = [leu]/ [leu] c and 
n T = [trp]/[trp] c : 

Leu - Trp FBR expansion velocity: Vl(?il) = vz,T(ni,), (1) 
Leu FBR Trp" expansion velocity: Vt(^t) = ^t^X^t), (2) 

where we have defined the threshold function: 

rpf \ \ x for x < 1 , s 

T(a;) = { 1 for x > 1 (3) 

Note that the velocities are expressed as function of concentrations 72l and nx felt by the cells, which can 
differ from the external amino acid concentrations tt-el and tt-et with which the agar medium was prepared, 
see below. 

Amino acid uptake. Under steady state growth conditions, the uptake of amino acids that are not 
produced by the cell must equal the amount used in the production of cellular biomass. Thus, the amino 
acid uptake rates must be proportional to the growth rates: 

Leu uptake rate of Leu - Trp FBR : K L (n L ) = k L T(n L ) (4) 
Trp uptake rate of Lcu FBR Trp - : K T (n T ) = k T T(n T ) (5) 

where and k^ are the maximal leucine and tryptophan uptake rates, respectively. 

Amino acid secretion. Sine our yeast strains are feedback-resistant in the production of leucine or 
tryptophan, they cannot regulate production in response to the available amount of the end product. We 
therefore assume that yeast secrete amino acids at a constant rate, 

Leucine secretion rate of Leu FBR Trp - : -Rl(^l) = r^, (6) 
Tryptophan secretion rate of Leu - Trp FBR : i?x("T) = '"t, (7) 

see also Ref. [20]. 

Amino acid dynamics. The dynamics of leucine and tryptophan with concentrations [leu] and [trp] 
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can be described by 

^[lcu] = r h &r-K h ([leu])c L +D a A[leu], (8) 
d 

— [trp] = r T c L - A T ([trp]) c T + D a A [trp]. (9) 

The first and second term describe amino acid secretion and uptake, respectively, by Leu - Trp FBR cells (sur- 
face concentration cl) and Leu FBR Trp~ cells (surface concentration ct). The last term describes diffusion 
of amino acids with diffusion constant D a . 

Yeast colony growth resembles a chemostat. As described in more detail in Section S5, only cells in 
a small active layer near the colony boundary are actively growing. The growth dynamics within the active 
layer in the reference frame of the boundary is reminiscent of growth in a chemostat, in which nutrient are 
added at a constant rate that is balanced by outflux of waste (medium and cells) [21]. For the active layer 
of a colony, there is continuous influx of nutrients (e.g. glucose) due to diffusive flux towards the colony. 
Cells are removed 'behind' the active layer, because cells in the colony interior are not growing and therefore 
effectively 'gone'. When the colony grows at constant velocity, nutrient influx and cell 'outflux' are balanced, 
so that the number of actively growing cells remains constant . Thus, we can describe cell dynamics in terms 
of the fraction 

/ = (10) 

of Leu~ Trp FBR cells of the constant surface carrying capacity c = cl + ct- As explained in Section S5, 
chemostat-like growth also means that the Laplace operator A can be decomposed (in a frame co-moving 
with the frontier) into angular and radial diffusion in a specific way, e.g. for leucine 

A[leu] = [leu] - d ([leu] - [leu] E ). (11) 

Here, the first term describes angular diffusion along the front coordinate x, with dx = Rd(f>. The second 
term replaces radial diffusion by effectively describing the radial diffusive flux into or out of the colony as 
dilution with rate d into medium with concentration [1cu]e (the concentration in the agar medium far away 
from the colony, i.e. the concentrations with which the medium was prepared). With Equations (10), (11), 
the diffusion equations (8,9) become: 

^ = ^a a ^-d([leu]-[lcu] E )+r L c(l-/)-A^ L ([leu]) C / (12) 
= D a 9 ^-d([trp]-[trp] E )+r T c/-ii- T ([tr P ])c(l-/) (13) 

Non-dimensionalization. These equations can be written in a non-dimcnsionalized form for the 
rescaled amino acid concentrations = [leu]/ [leu] c and nx = [trp]/ [trp] c by expressing time in units 
of the inverse diffusive rate d and space in units of the diffusion length scale \J D a /d: 

K-n EL )+p L (l-/)-K L T(n L )/ (14) 



dt dx 2 



(n T - n ET ) + p T f - k t T(ti t )(1 - /) (15) 



dt dx 2 

with i =td and i = xj \J D a /d. The secretion and uptake parameters are now dimensionless: 

r L c r T c fc L c k T c 

PL = -Fj — T—j, Pt = FT ■) k l = r, — T—j, kt = ft i (16) 

[leu] c a [trpj c a [leuj c a [trp] c a 

Intuitively, in the reduced parameter c/([leu] c d), the 'free' secretion rate tl is reduced by the factor d 
due to diffusive loss of amino acids from the colony. In addition, it is rescaled from cellular to amino acid 
concentrations with the 'conversion factor' c/[leu] c . 



7 



Note that for our mutualistic yeast strains, the tryptophan secretion rate Pp, uptake rate /ct, and crossover 
concentration [trp] c are all lower than the corresponding quantities for leucine by a factor of 8, see Sec- 
tion SI. 4. Since all dimensionless rates depend on the ratio of a secretion or uptake rate and the crossover 
concentration, this means that the dimensionless secretion and uptake parameters have the same values for 
the two strains, 

rcscaled nutrient uptake rate: n = kl = ktj (17) 
rescaled nutrient secretion rate: p = Pl = Pt- (18) 

The numerical values of these dimensionless parameters are listed in Table S2. 



S2.2 Steady-state amino acid concentrations 

Steady-state amino acid concentrations. The cell concentrations cl and ct vary in angular direction 
due to spatial demixing into patches or sectors. If only one cell type is present, or on length scales large 
compared to the patch or sector width, the cells in the active layer are essentially well-mixed, and the 
problem becomes radially symmetric, i.e. independent of the radial angle <fi. In this case the amino acid 
dynamics of Equations (14) and (15) become 



1 d 
1 driT 



p(l - /) - reT(n L ) / - (n L - "el), 
pf- nT(n T ) (1 - /) - (n T - n ET ). 



(19) 
(20) 



For a calculation that takes into account cellular patches along the colony front, see Section S4. 

Since the nutrient dynamics are fast compared to colony growth, we can assume that the cell fraction / 
is constant on nutrient time scales. The nutrient dynamics equations Equations (14) and (15) then have the 
quasi-stationary solution 




"EL+P(1-/) 
1+Kf 

"EL + P(l -/)-«/ 

»ET+P/ 



™ET + Pf - «(1 - /) 



n m h(f) 
ndL(f) 

n mT (/) 

HdT(f) 



for / > fci 
for / < fci 

for / < /c2 
for / > fc2 



(case ?il < 1) 
(case ni, > 1) 

(case nx < 1) 
(case i%t > 1) 



(21) 
(22) 



with the limiting allele frequencies 

fc 



p-l+n EL 
1 P+k 
f = K+l-rt ET 



2» E L-(2 + K — p) 

2( P +k) 
2ti E t-(2+k-p) 
2(p+k) 



(23) 



In view of Equations (1) and (2), this means that the growth velocities of the two strains are given by: 



«L"mL(/) 



for / > fci 
for / < fci 



and 



Vr(f) 



vt n mT (/) 



for / < fc2 
for / > fc2 



(24) 



Crossover concentrations from linear to saturated growth. For the special case of single-strain 
colonics, consisting e.g. of only Leu - Trp FBR cells (/ = 1), the steady-state leucine concentration is 



"L 



"EL 
1+K 

"EL 



for TlEL < 1 + K, 

for n EL > 1 + k, 



(25) 



This means that the colony expansion velocity from Equation (24) as function of the external amino acid 
concentration is 



Vcolony = Vl(tIel) 



1 + K 



for n EL < 1 
for n E L > 1 



ft, 
re, 



(26) 



Thus, when measuring the colony velocity as function of the external amino acid concentration, as in Sec- 
tion SI. 3, the velocity passes from the linear regime (proportional to the amino acid concentration) to the 
saturated regime (plateau velocity) at uel = 1 + k. This threshold is higher than the crossover value of 
1 naively expected from Equation (1) because the amino acid concentration felt by the colony is lower 
than the external concentration by the factor (1 + k). Similar equations hold for single-strain colonies of 
Leu FBR Trp~ cells. In particular, the factor (1 + k) is the same for both strains, which we indeed observe 
experimentally (Fig. S2). 



S3 Generalized Moran model for mutualism at the front 

Allele frequency dynamics. The expansion of a circular yeast colony on an agar surface is a two- 
dimensional process. However, since only cells in a small active layer near the colony boundary are actively 
growing (Section S5) one can approximate the expansion by a one-dimensional process along the front 
coordinate x [22]. In addition, because colony front growth is similar to growth in a chemostat (Section S2.1), 
the effective population size remains approximately constant as long as the inflationary effect of the radius 
increase is not too large. In the context of population genetics, our two different yeast strains can be 
viewed as representing two different species or as two different alleles. In total, the selection dynamics 
for our two 'alleles' can be described within the framework of the Moran model, with the addition of a 
frequency-dependent growth rate and spatial diffusion [22-24]: 

2 



%= D ^ + <<0 /(I - /) + y/D R (f)f(l-f)T(x, r), (27) 

where f(x,r) is the fraction of Leu - Trp FBR cells as function of the coordinate x along the front and the 
time r measured in generations. The first term on the r.h.s. of Equation (27) describes spatial diffusion with 
diffusion constant D s , which arises because dividing yeast cells push neighboring cells around. The second 
term describes selection with the selection coefficient 

s(f) = VUfl-VrV) (2g) 

Selection arises due to differences in the growth velocities Vl and Vt of the two strains, as given by Equa- 
tion (24). Since the growth velocities depend on the 'allele frequency' / because of the mutualistic interaction, 
the selection coefficient s is frequency-dependent. The last term in Equation (27) describes genetic drift, i.e. 
noise due to cell number fluctuations. T = T(x,t) is an Ito delta-correlated Gaussian noise, and the genetic 
diffusion constant 

n (f) -n W) + W) m) 

characterizes the strength of the genetic drift. It is frequency-dependent due to its dependency on the growth 
velocities. However, its frequency-dependence is weak compared to the frequency-dependence of s(f) so that 
we use D g (f) « D g . 



S3.1 The effect of mutualistic selection 

We first examine the effect of the mutualistic interaction which is incorporated in the selection coefficient 
s(f). If the selection term dominates in Equation (27), for example in a well-mixed culture with high cell 
numbers (so D g « 0), the dynamics become 

§£ = *(/)/(!-/)■ (30) 

As shown above, amino acid secretion and uptake is symmetric for our two strains because of k — kl = kt 
and p = pL = Pt (Equations (17) and (18)). For simplicity, we neglect here the 2% difference in growth 
velocities and use v = i>l = i>t- (We will consider the effect of the small growth velocity difference later in 
Section S3. 2.) With these assumptions, the selection coefficient can be obtained from the steady-state amino 
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acid concentrations (21,22) via the growth velocities (24). Its form depends on which of the limiting allele 
frequencies fci and fc2 given in Equation (23) is larger, 



fci < fc2 & "el + n ET < 2 
In the case fci < /c2, the selection coefficient is: 



«(/) 



/ + (l-/)"mT(/) 

" mL (/)-» m T(/) 

/n mL (/)+(l-/)n mT (/) 

"m L (/)-l 

/n m L(/) + (l-/)n mT 

and in the opposite case fci < /c2 : 

1-".t(/) 
/+(l-/)n mT (/) 

«(/) = <! o 

"mL(/)-l 

/« mL (/) + (l-/)r lmT 



for 0 < / < f cl 
for /ci < / < /C2 
for /C2 < / < 1 



for 0 < / < f C 2 
for /c2 < / < /ci 
for /ci < / < 1 



(regime Ct with n E > 1,Mt < 1) 
(regime M with n E < 1,71t < 1) 
(regime Cl with n E < l,n,T > 1) 



(regime Ct with til > 1,^t < 1) 
(regime N with ?t,l > 1, Ht > 1) 
(regime Cl with til < 1,tit > 1) 



(31) 



(32) 



(33) 



In these equations, we have labeled the frequency region in which both strains are amino-acid limited 
(m,, rix < 1) as M for 'mutualistic', since in this regime both strains benefit from the additional amino 
acids due to the mutualistic interaction. This contrasts with the 'neutral' N region in which both strains 
are saturated for their required amino acids, n El rt > 1, so that there is no benefit due to mutualism. The 
'commensal' regions Ct and Cl denote the cases (nL > 1, tit < 1) and (til < 1, tit > 1), respectively, in 
which one strain is amino-acid limited but the other is not. Some of these regions may vanish depending on 
whether the limiting fractions fci and fc2 are in the relevant interval (0, 1): 



0 < fci < 1 

1 > h2 > 0 



<(=> 1 - p < n EL < 1 

1 — p < ?lET < 1 



(34) 
(35) 



'Phase diagram'. In total, the selection coefficient s(f) depends on the parameters K, p, and the 
external amino acid concentrations n EE , tiet- Depending on the values of the allele frequencies fci and fc2 
(that mark the transitions between the commensal and the mutualistic or neutral regimes of the selection 
coefficient) in Equations (31,34,35), s(f) could be in one of ten possible 'phases', M, CtM, MCl, CtMCl, 
CtNCl, CtN, NCl, Ct, Cl, and N. Here, for example, the notation CtM means that s(f) has the 
functional form of the Ct regime (tryptophan limiting, leucine not) in Equation (31) for 0 < / < fci and 
the functional form of the M regime (tryptophan and leucine both limiting) in Equation (31) for fci < / < 1. 
As we will see below, for yeast parameters, only five of the ten possible phases appear; they are shown in 
the (helmet) plane in Fig. S3A. To interpret the graphs of the selection term s(f) /(l — /) as function of 
the allele frequency /, note that positive values of the selection term mean that selection acts to increase /, 
while negative values drive towards smaller /. A selection term of zero, s(f) f(l — f) = 0, defines a fixed 
point of the selection dynamics. The fixed point is stable if the selection term is positive to the left and 
negative to the right of it. Thus, in the CtMCl phase, selection drives the system towards a stable fixed 
point within the M regime (unless / = 0 or / = 1, which are absorbing boundaries of the system). In the 
CtNCl, CtN and NCl phases, selection drives the system from the commensal regimes Ct or Cl towards 
the neutral N regime where the selective force is zero. If genetic drift is added, it will dominate in this 
neutral regime, causing random fluctuations of the allele frequency /. 

The phase diagram in Fig. S3A shows that for all amino acid concentrations below the crossover con- 
centration, tt-el, Het < 1 + k, selection acts towards the interior of the frequency interval (0,1) and thus 
promotes mixing (red area in Fig. S3(A)). For very small amino acid concentrations (dark red triangle), 
selection drives the system towards a well-defined stable frequency /* , while for intermediate concentrations 
(light red area), selection only drives towards a neutral region within the interval (0, 1). For larger concen- 
trations, nEL, fiET > 1 + there is no barrier to fixation at / = 0 or / = 1. If genetic drift is added, random 
fluctuations in / will lead to fixation at one of these absorbing boundaries. The change in the selection term 
s (/)/(l — /) along the 'diagonal' n E = n E L = "et is shown in Fig. S3B. 

Mutualism with a fixed point. We first take a closer look at cases in which s(f) has a mutualistic M 
region, where both amino acid concentrations are below the cross-over concentration required for maximum 



10 




Figure S3: 'Phases' of mutualistic selection for yeast parameters as listed in Table S2. A) Deterministic 
'phase diagram' in the plane of external amino acid concentrations (tiel, J1et)- The insets show the selection 
term s(f) /(l — /) in Equation (27) on a vertical scale of —0.5 to 0.5 as a function of the fraction / of 
leucine-requiring cells in the interval 0 < / < 1. A positive selection term, s(f) f(l — f) > 0, acts to increase 
/, while negative values drive towards smaller /. The selective force vanishes for a selection term of zero, 
s (f) /(I — /) = 0- The selection term can be in one of five possible 'phases' designated as CtMCl, CtNCl, 
CtN, NCl, and N. In the phase names, M denotes the mutualistic regime in which both strains are amino- 
acid- limited, til, tit < 1, so that they benefit from their mutualistic interaction. N represents the neutral 
regime with non-limiting amino acid concentrations til, tit > 1 and therefore no mutualism. Ct denotes a 
commensal regime in which the leucine requiring strain is saturated, til > 1, while the tryptophan-rcquiring 
strain is still limited, tit < 0. Cl describes the inverse situation with Ul < 1, tit > 1- The vertical red 
and green lines in the insets mark the transition frequencies fci and fc2 between the different regimes. The 
red shades phases CtMCl and CtNCl represent mutualistic phases, in which selection drives the system 
towards a stable frequency /* £ (0, 1) within the M regime (dark red) or towards a neutral N region in the 
interior of the interval / £ (0, 1) (light red). In the green shaded, non-mutualistic phases CtN, N, and NCl, 
the system has no barrier to fixation at / = 0 and/or / = 1. B) Selection term s(/)/(l — /) along the 
diagonal ttel = ttet = tie- The vertical red and green lines in the insets mark the transition frequencies /ci 
and fc2 as in (A). The black horizontal lines at ±s c , with s c = D^/D s = 0.18, become comparable to the 
selection term around tie = 0.25. 



growth, i.e. tt-el + ttet < 2 + k — p, and tjeLj tiet < 1 + «. This mutualistic regime is characterized by a 
uniquely determined fixed point /* for which s(/*) = 0: 

f* = \ + \ o T \l + ( 2 T T ( n EL - tiet) (36) 

2 2 2p(l + k) + K{n EL + tiet) 

Because of the symmetry of the mutualistic interaction parameters (k = kl = «t, P = Pl = Pt), the station- 
ary allele frequency /* equals 1/2 unless driven away from it by unequal external nutrient concentrations 
tiel 7^ tiet- Vanishing selection s(f) = 0 in the mutualistic regime means that the amino acid concentrations 
til(/*) = tit(/*) are equal, see Equation (32): 

n L (f ) = n T (f ) = . (37) 

As expected, the amino acid concentrations increase with the external concentrations tiel, tiet, and the 
cellular secretion rate p, but decrease with the uptake rate k. 

Obligate mutualism. In the obligate case, ttel = tiet = 0, the steady-state allele frequency is 
/* = 1/2 due to the symmetry, as observed experimentally. The steady-state amino acid concentrations 
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Figure S4: A) The frequency-dependent selection term s(/)/(l — /) of Equation (27) for different amino 
acid asymmetries ?iel ~ ?iet , with fixed total amino acid amount tiel 4- ?^et = 0.6. For increasing asymmetry 
n-EL — n-ET, fixation of leucine consumers at / = 1 becomes more likely, because (i) the stable fraction /* 
(indicated by vertical lines) of leucine consumers moves closer to / — 1, and (ii) the barrier towards fixation 
at / = 1 (the absolute value of the minimum of the selection term) decreases. B,C) Dynamics of the colony 
boundary fraction /, as function of radial distance R — Rq from the homeland of radius Rq, for colonies with 
different inoculation fractions fo- B) For obligate mutualism, the the boundary fraction reaches the steady- 
state fraction /* = 0.5 within a 0.5mm radial expansion from the homeland with weak transient oscillations. 
C) For no mutualism, the boundary fraction increases during expansion because blue Leu~ Trp FBR cells are 
2% more fit than yellow Lcu FBR Trp - cells. Data points in (B,C) are from flow cytometry of cells harvested 
from colony boundaries. Solid lines are solutions to selection dynamics with selection coefficient s(f) of 
Equation (28) for obligate and s = 0.02 for no mutualism. 



are n^(f*) = Ht(/*) = p/(2 + k), so that the expansion velocity according to Equation (24) is vn^,(f*) = 
vtit(/*) = v p/(2 + k). In our experiment, obligate mutualists expand with a velocity v/2 that is half the 
maximum velocity, sec Figure 4F in the main text. Thus, we have 

P = 1+ f • (38) 

We use this condition to determine the value of the rescaled secretion rate p from the value of k. Equation (38) 
shows that for our mutualistic strains p > 1, so that the conditions n^L, "et > 1 — p in Equations (34) and 
(35) are always satisfied. In consequence, only the phases CtMCl, CtNCl, CtN, NCl, and N are relevant 
for our system, as shown in Fig. S3. 

Absorbing boundaries for obligate mutualism. In the obligate case tiel = ^et = 0, the selection 
coefficient in the mutualistic M regime in Equation (32), 

diverges when the fraction / of leucine-requiring cells approaches 0 or 1. The reason is, for example in the 
case / — > 0, that the growth velocity Vl of the rare Leu" Trp FBR cells approaches a constant, while the 
velocity Vr of the abundant Leu FBR Trp _ cells vanishes linearly with /. Similarly, the average population 
velocity /Vl + (1 — f)Vr approaches zero linearly with /. Thus, as / —> 0, population growth ceases and 
thereby prevents the system from reaching the boundary / = 0 within finite time. A population with / > 0 
will therefore never fix at / = 0. Note that the selection term s(/)/(l — /) in Equation (27) remains finite 
for all / e [0,1]. 

Asymmetric mutualism. If e.g. more leucine than tryptophan is supplied in the medium, uel — tiet > 
0, mutualism is asymmetric, favoring leucine consumers, see Fig. S4. Compared to symmetric mutualism 
with a stable fraction /* = 0.5, the stable fraction is closer to / = 1, and the barrier towards fixation at 
/ = 1 (the maximum value of the absolute selection term in the interval [/*, 1]) is smaller. Both effects make 
fixation of leucine consumers at / = 1 more likely. Likewise, tryptophan consumers are more likely to fix if 
there is more tryptophan than leucine in the medium. 

Mutualism with a neutral region. The phase diagram in Fig. S3 exhibits an unusual phase CtNCl 
in which the selection coefficient s(f) = 0 for a finite region in the interior of the interval (0, 1). In this 
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phase, selection in the commensal Ct and Cl regimes drives the system towards the neutral N regime, but 
dynamics within this regime is neutral. The existence of the neutral N regime is due to the plateaus in 
the velocity curves in Fig. S2: If both amino acid concentrations are Ul, «t > 1, the precise concentration 
values do not matter since the colony always grows at maximum velocity. 

No mutualism. For uel, ^et > 1 + «, the system does not have a barrier towards fixation at one or 
both of the absorbing frequency boundaries / = 0, =1 (green regions in Fig. S3). In the NCl and CtN 
phases, selection acts towards the neutral N regime in which s = 0. Since the neutral regions connect to 
/ = 0 and/or / = 1, number fluctuations can push the system towards an absorbing boundary, thus leading 
to demixing. In the N phase, genetic drift is the only acting force, leading to eventual fixation at / = 0 or 
/ = !■ 



S3. 2 The effect of small fitness differences 

So far we have neglected the small fitness advantage = 0.02 of the Leu~ Trp FBR strain over the 
Leu FBR Trp _ strain (Section SI. 2). For low external amino acid concentrations, mutualism is a strong 
selective force with selection coefficients large compared to s&, see Fig. S3, and it is justified to neglect this 
fitness difference. Experimentally, this approximation results in the observed symmetric steady-state cell 
fraction /* = 0.5 'despite' the fitness difference. However, for higher amino acid concentrations, the fitness 
difference becomes relevant as mutualistic selection becomes smaller. 

No mutualism. For no mutualism, the fitness difference is the only selective force so that the dynamics 
of the fraction / of leucine-requiring cells becomes 

^ = Sd /(l-/). (40) 

Thus, the frequency changes from the start frequency fo as 

fn e SdT 

W = i + t(s dT n ( 41 ) 

1 + /o(e SdT - 1) 

and approches fixation of Lcu~ Trp , i.e. / = 1, for long times. This equation has been used to predict 
the cell fraction dynamics in Figure 3B in the main text and Fig. S4C. 

Weak mutualism. The fitness difference is also important for weak mutualism, in particular when the 
selection coefficient has neutral regimes, as in the N, CtN, NCl, and CtNCl phases: It 'tilts' flat the neutral 
regimes towards / = 1. This results in selection towards fc2 in the CtNCl phase, increases the likelihood 
for fixation at / = 1 in the N and CtN phases, while it decreases the probability of fixation at / = 0 in the 
NCl phase. This bias has been taken into account in the theoretical phase diagrams shown in Figure 4B in 
the main text, and explains at least part of the asymmetries in the experimental phase diagram in Figure 
4A,B in the main text. 



S3. 3 Transient approach to steady-state growth 

Fig. S4B,C show the dynamics of the boundary fraction / of blue cells as function of the radial distance from 
the homeland, which in our spatial expansion is analogous to time. In the case of obligate mutualism, the 
experimental boundary fraction converges to the steady-state value /* = 0.5, independent of the inoculation 
fraction /q. The approach happens rapidly, within less than 1mm distance from the homeland with a small 
overshoot towards higher /, a feature expected for cross-feeding interactions [25]. Our simple mutualism 
model, designed for steady-state colony expansion, fails to produce this overshoot (although it does capture 
the time scale of the approach to steady-state, see solid lines in (B)), as well as asymmetric time lags to 
initiate growth (colonies with fo = 0.15 are smaller since they take longer to start growing, see Section SI. 3, 
and colonies with fo = 0.01 do not grow at all). 

In contrast to obligate mutualism, selection for no mutualism is not frequency-dependent, and the colony 
boundary fraction / depends on the start fraction fo (Fig. S4C). The fraction / of blue cells increases 
during expansion because the blue tryptophan producers have a 2% fitness advantage over the yellow leucine 
producers under these conditions (Section SI. 2). 
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S3. 4 Antagonism between selection and genetic drift 



As shown in Section S3.1, the mutualistic selection coefficient selects for 'mixing', i.e. for cell fractions 
/ in the interior of the interval (0,1), as long as the amino acid concentrations are below the crossover 
concentrations Bel, iet = 1 + as = 1.18. However, experimentally local demixing into domains with / = 0 
or / = 1 occurs for much lower concentrations jt-elj^et > 0.25, see Figure 4 in the main text. This is due 
to the genetic drift term in Equation (27). In order to compare the strengths of mutualism and genetic 
drift, we non-dimensionalized Equation (27) by measuring distance in units of D s /D g and time in units of 
s c = D s /Dl to obtain [22,24] 

S = 0 + 4r /(1 - /)+ ^ r ^ r(J ' ?) - (42) 

with the dimensionless distance x = xj (D s /D g ) and time r = tj (D s /D^) As shown in Ref. [24], the strength 
of mutualism and genetic drift become comparable if mutualistic selection becomes comparable to the critical 
selection coefficient 




(43) 



that describes the strength of local demixing due to genetic drift. The 'mutualistic barrier' to demixing, 
i.e. fixation at / = 0 or / = 1 when the system is close to the stable fraction /*, can be estimated as the 
maximum of the selection term |s(/)/(l — /)|, see Fig. S3B. One can see in Fig. S3B and in Figure 4D in 
the main text, that for symmetric mutualism with tie = Hel = ^et this barrier becomes comparable to s c 
at about tie = 0.25. Thus, we expect that mutualism is dominant for ue < 0.25, leading to expansion in a 
mixed pattern, while genetic drift dominates for > 0.25, leading to demixed expansion. This is indeed 
found experimentally, see Figure 4 in the main text. 

For asymmetric mutualism, local fixation of the favored strain becomes more likely because the stable 
cell fraction /* moves closer to a fixation boundary / = 0 or / = 1, and because the 'mutualistic barrier' to 
fixation at that boundary becomes smaller, see Fig. S4. A more sophisticated analysis of the critical selection 
coefficient s c for asymmetric mutualism is the topic of a forthcoming theoretical paper [26] . 



S4 Model for patches of obligate mutualists 

In colonies of obligate mutualists, the two mutualistic strains form a pattern of blue and yellow 'patches', 
see Figure 2 in the main text. Because of the effective symmetry of our mutualistic interaction, yellow and 
blue patches have the same width, see Fig. S5A. After an initial transient of 1-4 days, the patch pattern 
becomes stationary, i.e. although patches may form and disappear, the overall pattern and patch widths 
remain statistically the same during the expansion. 

In this section, we derive an approximation for the characteristic width of the patches during steady-state 
growth by considering the amino acid dynamics. Basically, the width of a patch of leucine-consuming cells 
is limited by how far leucine can diffuse into it from a neighboring leucine-producing patch before cellular 
uptake makes its concentration too low for cells in the patch interior to grow well. Although the patches are 
probably not completely demixed, here we first assume, for simplicity, that a yellow patch consists entirely 
of yellow Leu FBR Trp~ cells, and a blue patch consists entirely of blue Leu~ Trp FBR cells. We denote the 
average width of a blue and yellow patches parallel to the front as Ll an d L^, respectively, see Fig. S5B. 

Diffusion with sources and sinks. We first consider leucine, whose dynamics in the active layer of a 
colony is described by Equation (14). In a leucine 'production patch' of Leu FBR Trp~ cells, we have / = 0 
and thus 

Id D a d 2 

Production patch: ddi UL=P ^ nh ~ " EL ^ + ~d 'dx 2 ™ L ' 

Similarly, the leucine dynamics in a leucine 'consumption patch' with / = 1 is given by 

Id £) a Q 2 

Consumption patch: ^^" L = ~ Kn L - («l - "el) + ~^f~g^2 nL - ( 45 ) 
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Figure S5: Determination of patch widths along the front coordinate x. A) Average patch widths (parallel 
to the front) as function of the radial distance R — Ro from the homeland with radius Rq, determined via 
different methods for the same colony of obligate mutualists. Yellow patch width (yellow data) and blue 
patch width (blue data) are the average angular local maximum-to-minimum and minimum-to-maximum 
distances, respectively, of the derivative of the yellow fluorescence intensity. The width of all patches (green 
data) is reproduced from Figure 2D in the main text. Since it is is calculated as the circumference divided 
by number of local maxima of the yellow fluorescence intensity, it docs not equal the average of the yellow 
and blue patch width. All patch widths are identical within error bars and saturate at L = 52 [im (black 
horizontal line) for radial distances larger than 1mm. B) The amino acid leucine is secreted by Leu FBR Trp - 
in yellow production patches of width Lt, and consumed by Leu - Trp FBR cells in blue consumption patches 
of size Ll- The leucine concentration til thus 'piles up' in the production patches, and decreases in the 
consumption patches. The width of the consumption patch is limited by this decrease, as explained in the 
text. 



Here, the uptake rate in the leucine-limited regime (n^ < 1) has been used; for higher amino acid concen- 
trations, the sink term would equal — kl- 

As shown in Fig. S5B, a consumption patch is surrounded on both sides by production patches. Since 
we assume all consumption and production patches to have the same widths Ll and Lt, respectively, the 
symmetry of the problem dictates that, in steady state, there is no amino acid current through the middle 
of the patches. Using the coordinate system shown in Fig. S5B, this leads to the boundary conditions 



d_ 

dx 



0 



and 



d_ 

dx 



»L 



0. 



(46) 



In addition, the amino acid concentration and current have to be continuous at the patch boundaries: 
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ox 



continuous at x = ±- 



(47) 



Solution to the diffusion problem. Since we are interested in describing the pattern of the mutualistic 
strains during steady-state growth, we solve the diffusion equations (44) and (45) for the stationary state in 
which the amino acid profiles do not change over time. The diffusion problem (44,45) with the boundary 
conditions (46,47) can be solved analytically in terms of hyperbolic functions. The minimal and maximal 
leucine concentrations within the consumption patch are attained in the middle x = 0 and at the boundary 
x = Ll/2 of the consumption patch, respectively. For obligate mutualism with tiel = 0, these values can be 
obtained from the full analytic solution: 



«L,min = n h (0) = 
max = n L (L L /2) = 



sinh 



sinh 



2^ 



cosh 



2L 



(48) 
(49) 
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where 



*-Vf a " d ! - = V^ (50) 

are the length scales characterizing the 'piling up' of leucine in the production patch and the decay of the 
leucine concentration in the consumption patch, respectively. We also used the abbreviation 



Velocity gradients within consumption patches. In the amino acid-limited regime, the growth 
velocity is proportional to the concentration of the required amino acid, see Equations (1,2). This means 
that the growth velocity varies within a consumption patch: It is smallest in the center and largest at the 
boundary. If cells grew only in the expansion direction, these velocity differences along the colony front 
would lead to non-uniform front-propagation and thus to an unstable, undulated colony front. However, 
experiments show that the obligate mutualists expand together with a smooth colony front, see Figure 2 in 
the main text. In addition, the velocity differences generate a selection force perpendicular to the expansion 
direction. In a leucine consumption patch, this perpendicular velocity gradient can be estimated by [6, 27] 



^J-L = vl yj nl >max - nl min . (52) 

There must be a force that counteracts the selective force generated by the velocity differences within a path. 

Patch boundary diffusion. Colony growth of non-motile yeast cells proceeds due to cell division, 
which results in pushing and shuffling of cells. This leads to a diffusion-like process that occurs parallel to 
as well as perpendicular to the expansion direction. Experimentally, this diffusive mixing can be observed as 
the diffusion of sector boundaries with diffusion constant D s /T g , where r g is the generation time at the patch 
boundary [27,28]. Since the growth velocity at the boundary equals v = i'L"-L,max, and since the growth rate 
at the boundary is g = 1/t s = v/b (Section SI. 3), we can estimate the diffusion 'velocity' required to diffuse 
on a length scale of half a patch of width as 

_ 4{D s /r g ) _ AD 2L 

v ^ = ^- = ^zi — = — u~ ' (53) 

where we have introduced the boundary diffusion length scale l\, = 2D s /b. 

Balance condition determines patch width. This diffusive process can counteract the imbalances 
created by the velocity differences if both velocities are of the same order of magnitude, 

V ±L = y DL - (54) 

An analogous equation holds for the width Lx of tryptophan consumption patches. In the case of obligate 
mutualism, i.e. no external amino acids tt-el = ^et = 0, the mutualistic interaction is fully symmetric so 
that Ll = Lt = L. In this case, Equation (54) yields 

L = 2J b cothfij , (55) 

which can be solved approximately for small patch widths to yield Eq. [5] in the main text, 

L = 2 v/^A. (56) 



S5 Radial flux of amino acids into and out of a colony 

In this section, we derive an approximate description for the diffusive flux of amino acids into and out of 
the actively growing boundary layer of a colony. We will show that, on the experimental timescales of many 
hours to days, the amino acid diffusion dynamics can be approximately described by a chemostat-like process 
as used in Eq. [1] in the main text, and we will derive an expression for the effective influx / outflux rate d 
in terms of the colony growth parameters. 
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Figure S6: A) Geometry to calculate amino acid diffusion within and out of a colony. In a colony of 
radius r c , only cells in an active layer of size b = r c — r\ near the colony boundary grow, while cells 
in the colony interior r < n do not grow due to depletion of essential nutrients (e.g. glucose). In the 
active layer, amino acids diffuse, and are produced and taken up by cells. Outside the colony, amino 
acids diffuse freely and create a diffusion zone of radial size ?'d — r c ~ \/AD a t. B) Radial amino acid 
concentration profile n^(r) = [leu] r /[lcu] c for yeast parameters (Table S2) for external nutrient concentrations 
™el = [leu]/ [leu] c = 0 (red) and n EL = 0.5 (blue). 

55.1 Nutrient diffusion and colony growth 

Diffusion time scales. The cells in a growing yeast colony obtain their nutrients, e.g. glucose and amino 
acids, from the agar medium on which they grow. Nutrients reach cells via diffusion with a diffusion constant 
D nn t of a few mm 2 /h [15,29]. . This parameter value means that diffusion can sample the 5 mm height 
of the agar in the Petri dish within a few generations (generation time r g « 1.5h). Thus, on experimental 
timcscales of several days, we integrate out the diffusion perpendicular to the agar surface and only consider 
the two-dimensional diffusion of nutrients in the plane of the colony, see Fig. S6. In addition, diffusion is fast 
compared to the colony expansion velocity of v sa 20/^m/h: a colony will outrun nutrient diffusion only after 
about D nut /v 2 w 1 year. We can therefore ignore colony expansion for nutrient dynamics on experimental 
timescales of several days. 

The active layer. We first consider the diffusion and uptake of an essential nutrient that is not produced 
by the cells, such as glucose. After an initial transient of 1-2 days, in which the colony depletes the glucose 
beneath it, new glucose reaches the colony by radial diffusion towards the colony. This glucose allows only 
the cells in the active layer of size b = 40 fjja near the colony boundary to grow [11, 12]. Consistent with the 
chcmostat-like nature of the frontier region, these cells in the colony boundary consume all the glucose, so 
that cells in the colony interior cannot grow. This geometry is shown in Fig. S6A. 

55.2 Radial amino acid dynamics 

We now consider the radial diffusion dynamics of leucine (tryptophan works analogously) , which is produced 
or taken up by cells in the active layer of a colony consisting of one or both of the mutualistic strains. 
Here, we are interested in length scales larger than the angular width of patches or sectors. We therefore 
consider the cells in the active layer to be essentially well-mixed, so that they can be described by the 
r-independent concentrations cl and ct of Leu - Trp FBR and Lcu FBR Trp - cells, respectively. The two- 
dimensional dynamics of the radial leucine concentration [lcu] r = [lcu](r, t) in the active layer can then be 
described by the radially symmetric diffusion equation, 



for n < r < r c : 



d_ 



[leu] r = D a - — 
r or 




[leu] c 



c L [leu] r + r L c T . 



(57) 
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where r c is the colony radius and r c — r\ = b is the size of the active layer. Outside of the colony, leucine 
diffuses freely, 

19,. , 19 / <9 \ . 

fOTr - rc: D- a di [leu]r = r^{ r d-r) [leu]r - (58) 

This diffusion causes a diffusion zone around the colony, within which the leucine concentration differs 
from its concentration n E L m the agar medium far away from the colony (the concentration with which the 
medium was prepared). On the experimentally relevant times of 4-7 days, the size of the diffusion layer, 
'"d — r c ~ y4JDaf, increases very slowly from 34mm to 44mm. For distances larger than this, the leucine 
concentration equals its value far away from the colony, 

[leu] r (r > r d ) = [leu] E . (59) 

In a quasi-stationary approximation, we can solve the diffusion equations (57,58) for their stationary states 
with the boundary condition (59). In addition to the boundary condition (59), the amino acid concentration 
and current have to be continuous at the connection between the active layer and the diffusion zone: 

. . 9[leu] r . . . 

leu L, — continuous at r — r c . (60) 

or 

Furthermore, since cells in the colony interior are inactive, amino acids diffuse freely for r < n, with a 
zero-current condition in the colony center at r = 0 because of symmetry. This translates into a zero-current 
condition at the boundary between the active layer and the colony interior: 

^(r = n) = 0. (61) 

Solution to the diffusion problem. The diffusion problem (57,58) with the boundary conditions (59- 
61) can be solved analytically. The concentration within the active layer is a superposition of the modified 
Besscl functions of the first and second kind, while the concentration decays logarithmically in the diffusion 
zone. It is plotted for yeast parameters in Fig. S6B. 

Flux into and out of the active layer. As can be seen in Fig. S6B, radial diffusion happens on length 
scales large compared to the radial size of the active layer b = 40 fim. Thus, we can approximate the leucine 
concentration within the active layer, n^, by the concentration at the colony boundary, [leu] = [leu] r (r c ). 
The flux at the colony boundary can be calculated from the exact solution to the diffusion problem and has 
the simple form 

([leu] -[leu] E ). (62) 



. _ . <9[leu] r 
JL = Jh(r c ) = D a 



dr 



.In 



Due to the boundary flux j'l , the number of leucine molecules in an area element b Ax of the active layer 
changes in a time increment At according to 

[[leu] (t + At) - [leu] (t)]bAx = j h Ax At, (63) 



or 



[leu] = 3 ± = -d ([leu] - [leu] E ) (64) 



dt l 1 b 
with the diffusive flux constant 

d = ^ r- = ^ - rr. (65) 

6r c ln(^j 6r c ln(l + ^j 

This diffusive flux is balanced by leucine secretion and uptake in the active layer, leading to an effective 
dynamic equation for the nutrient concentration Hl in the active layer 

d 

— [leu] = -fc L c L [leu] + r L c T - d ([leu] - [leu] E )- (66) 

A similar equation holds for tryptophan dynamics. With the amino acid diffusion constant D a = 3 mm /h, 
the size b = 40 /im of the active layer, a colony radius of r c = 3 — 5 mm, and the size y/AD^t ~ 34 — 44 mm 
of the diffusion zone, the diffusive flux constant is of the order of d ~ 5/h. 
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